clc;
clear;
f=@(x)exp(x.^2).*(x<=2)+...
    80./(4-sin(16*pi*x)).*(x>2);

[x1,f1]=intfunc(f,0,4,100);

plot(x1,f1,x1(end),f1(end),'o'),
% x1(end), f1(end), 
% 'o' 用于在图上标记最后一个采样点 (x1(end), f1(end))，
% 即在最后一个采样点处绘制一个圆圈。
f1(end)%f1(end) 用于获取数组 f1 的最后一个元素，即最后一个采样点处的函数值。